Sum-rules for electron energy-loss near-edge spectra 
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We derive four sum-rule expressions for spectra measured in electron energy-loss near edge struc- 
ture experiments. These sum-rules permit the determination spin and orbital magnetic moments, 
spin-orbit interaction and number of states, analogously to the sum rules of x-ray magnetic circular 
dichroism. The derivation of the sum-rules is based on dynamical electron diffraction theory and the 
properties of the mixed dynamic form-factor. The accuracy of the sum-rules is tested by a complete 
evaluation of the thickness dependent electron energy-loss spectra for iron, cobalt, and nickel crys- 
tals. We find that the sum-rules reproduce both spin and orbital moments with very good accuracy. 
Our results provide a foundation for the use of the energy loss magnetic chiral dichroism technique 
as a quantitative probe of element specific magnetic properties. 

PACS numbers; 68.37.Lp, 75.75.-|-a, 75.40.Mg 
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X-ray magnetic circular dichroismiii^'^ (XMCD) is a 
powerful and widely used spectroscopic technique for 
the investigation of magnetic materials. The strength 
of the XMCD stems particularly from the well-known 
sum-rules, which enable one to obtain element specific 
information on spin and orbital magnetic moments in 
the material from an integration of the measured XMCD 
spectrum4ii^i^ 

The existence of the electron analogue to the X-ray 
radiation based XMCD was recently demonstrated by 
Schattschneider et al^ An electron beam passing through 
a magnetic material in an electron transmission micro- 
scope (TEM) was proved to exhibit a magnetic dichroism 
phenomenon, which has been named energy-loss mag- 
netic chiral dichroism (EMCD). The EMCD technique 
has the advantage that, due to the short de Broglic 
wavelength of the electrons, an excellent lateral resolu- 
tion can be achieved. Modern TEMs easily reach sub- 
nanometric lateral resolutions, though this has not yet 
been achieved for EMCD. By improvements of the ex- 
perimental geometries^i^ the signal-to-noise ratio (SNR) 
of the EMCD has recently been significantly enhanced. 
Lateral resolutions of 10 nm and below have been reached 
with a very good SNR, exceeding thereby the resolution 
accomplished in current XMCD experiments. 

In spite of the encouraging recent achievements of the 
EMCD, there is a lack of quantitative interpretation of 
the measured spectra connecting these to magnetic prop- 
erties of the studied material. Here we fill this gap by de- 
riving spin and orbital momentum sum-rule expressions 
for the EMCD technique. Our results thus provide a for- 
malism to establish this new experimental method as a 
quantitative probe which provides element specific infor- 
mation on spin and orbital magnetism. 

EMCD is essentially an electron energy-loss near-edge 
structure (ELNES) experiment, in which a core state 
electron is excited to an unpopulated valence state above 
the Fermi level. Precise evaluation of the double dif- 
ferential scattering cross-section [DDSCS, <j{E)], which 



is an output of ELNES experiment, requires the combi- 
nation of dynamical diffraction theory and microscopic 
calculations 

In the following we present the derivation of the EMCD 
sum-rules. It was demonstrated that the dichroic effect is 
related to the imaginary part of the mixed dynamic form 
factor [MDFF; usually denoted as ^(q, q', E) - a. function 
of energy and momentum transfer vectors] f^i^ Because of 
its hcrmiticity, the imaginary part of the MDFFs can be 
separated out by performing two measurements - simi- 
larly to XMCD - in which the order of q, q' is reversed. 
It is assumed that both measurements are performed in 
symmetrical geometry, so that the dynamical diffraction 
equations lead to the same set of Bloch coefficients and 
Bloch-wave vectors. For example, a reversal of the sam- 
ple magnetization without changing the experimental ge- 
ometry would provide such two measuremcntSFi^ Then 
the difference of two such spectra, let's denote them by 
CT+ (E) and ct_ (E) , is determined solely by the imaginary 
parts of contributing MDFFs. The sum of these spectra 
is a function of the real parts of MDFFs. 

The MDFF itself is given by the following expression 
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where the sum is taken over initial (core) states and final 
(unoccupied) band states. The sum-rules are expressed 
in terms of integrals of the signal from a particular edge 
over the whole energy spectrum. Owing to the ^-function 
the integration of the MDFF over energy is trivial and 
the sum over final states can be formally expressed as 
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where p is the ground state density matrix. 

The next step is to consider only dipole allowed allowed 
transitions I ^ I + 1 = L (e.g. 2p —>■ 3d in transition met- 
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als). All such transitions are contained in the second 
term of the Rayleigh expansion of the exponential e*'' '", 
which is equal to 3ji(qr)q • r/{qr) [j\{x) is a spherical 
Bessel function of order A] . For small qr this expression 
reduces to a simple scalar product q • r appearing also in 
the Taylor expansion of the exponential (and often used 
in XMCD interpretation). But for qr > 1, which is a 
common situation in ELNES, these expressions arc very 
differentii and the Rayleigh expansion provides signifi- 
cantly higher accuracy. The second assumption is that 
the radial part of the core state wavefunction is spin and 
j-independent. We note that an exact sum rule for the 
orbital momentum in the XMCD has been derivedji^ and 
in Refs.[H3it was demonstrated how to include spin and 
j-dependence. In principle these approaches can be ex- 
tended to the EMCD sum-rules derived below, but for 
clarity of presentation we will assume that these approx- 
imations are accurate enough. 



elements and summation over /i, /i' we obtain 

^Im[5(q,q')] = hlLiq, q'){2L - 1)Q ■ {I) (4) 



where 

ML{q,q') 
and Q = q X q'. 
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An EMCD analogy of the XMCD spin sum-rule is 
i Im[5(q,q')]|,,+ - Im[5(q, q')]|,,_ = 
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Under these assumptions the integrated MDFF (we use 
the same notation for energy-integrated MDFF omitting 
the energy parameter) becomes 

5(q,q')«127ri^0-i(g)),i0-i(g'));L 
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where I and L are the orbital quantum numbers of ini- 
tial and final states, respectively, pL is the L-projected 
part of the density matrix, {ji{q))iL is the integral of the 
spherical Bessel function, ji{qr), and radial parts of ini- 
tial and final wavcfunctions, and Yj^{q) are the spherical 
harmonics, q ~ q/q. The Wigner Sj-symbols appear due 
to integrals of spherical harmonics in the real space (for 
details see Ref . Hlh . Kronecker deltas S,i± select the par- 
ticular edge (J = I ± i). This expression is an ELNES 
analogy of Eq. (5) in Ref. @. The formal similarities of 
this expression and matrix elements for the x-ray absorp- 
tion spectra allow to use the same L-operator algebra as 
in Ref. d- Note, however, that in contrast to XMCD, 
where the particular photon polarization fj. can be cho- 
sen in a particular experiment, in the ELNES expression 
every MDFF term is a combination of all three "polariza- 
tions" including cross-terms between them (non-diagonal 
terms in fi ^ p,'). 

The orbital sum-rule is given as a sum of the integrated 
dichroic signal over both J = I ± ^ edges. In such case 
the last-but-one line of Eq. Q (the only J-dcpendent 
part) reduces to (2L — 1). After evaluation of all matrix 



where T is the the magnetic dipole operator. 

Equations (g]) and (O fully disclose the information 
contained in the energy integral of the imaginary part of 
MDFF, which can be re-formulated as 
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depending on particular edge J = ^ ± i . 



Because the real part of the MDFF does not change 
upon time inversion (in the dipole approximation, other- 
wise it only holds if the system also has inversion symme- 
try), it must be an even function of magnetic moments. 
It is possible to derive two other sum-rules for the real 
part of the MDFF. The first one is an analogy to the 
so called N sum-rule,— corresponding to the number of 
states 
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This is the only sum-rule, where a nonzero contribution 
arises from the unity operator in Eq. ([3]) - leading to the 
first term of Eq. ([5]). The second term is proportional 
to the trace of the density matrix, i.e. the number of 
i-clcctrons, Ne = {1} = Tt[pl]. 



The second sum-rule for the real part of the MDFF is 
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an analogy of the spin-orbit sum-rule^ 



and spin sum-rule 



i Rc[5(q,q')]|^+ - Rc[5(q,q')]|,,_ = 

- 2A/i(g,q') [(q.q')(L.S) + 

((q.L)(q-.S) + (q-.L)(q.S)) 
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where the dominant term is proportional to the mean 
value of the spin-orbit operator. 



In both these sum-rules there are some extra terms, 
which do not appear in their XMCD analogues. It is 
possible to show that these vanish in the limit of a 
strong enough magnetic field, where all magnetic mo- 
ments are aligned along the z-axis (magnetization axis) 
and small energy losses, when the z-component of mo- 
mentum transfer vectors vanishes. 



These four sum-rules, written in a rotationally invari- 
ant form, are a complete set of ELNES sum-rules within 
the dipole approximation. Ankudinov and Rehr— derived 
six X-ray absorption spectroscopy (XAS) sum-rules, be- 
cause there it is possible to prepare incoming light in a 
particular polarization — or 0). In the ELNES exper- 
iment every measurement is in general a combination of 
all possible polarizations including 'cross-terms'. Having 
two symmetric measurements, it is possible to set apart 
the signals which originate from the real or imaginary 
parts of the MDFF's. By separating the information ob- 
tained for the two edges we have four independent inte- 
grals of signals, which allow to construct the above-given 
four sum-rules. Expressed in simplifying terms, one can 
say that (i) the ELNES N sum-rule is a combination of 
XMCD N and anisotropic orbit sum-rules and (ii) the 
ELNES spin-orbit sum-rule is a combination of XMCD 
spin-orbit and anisotropic spin magnetic*^ sum-rules, re- 
spectively. 



To calculate the total DDSCS, we need to evaluate a 
sum over many MDFFsiii It is important to note that 
independent of the values of q, q', the MDFF is always 
given as a mean value of particular operators appearing 
in sum-rules. The summation then introduces only a 
prefactor (although a rather complex one) . The complete 
orbital sum-rule is 
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where Act = j[a+{E) — a-{E)]dE is an integrated dif- 
ference signal from a particular edge. All quantities en- 
tering the prefactor are accessible either from dynami- 
cal diffraction theory [Y^^J^^, is a product of Bloch co- 
efficients, Tjij'i' depends on sample thickness and Bloch 
wave- vectors] or from simple electronic structure calcula- 
tions [factor Ml{ci, q'), Eq. ^] ; they can be determined 
with good accuracy. More details can be found in our 
recent publication, Ref. [ill. 

If we assume that Qx,y <S Qzi which means that we 
neglect the z components of q, q' (physically this means, 
that the electron energy-loss is negligible when compared 
to acceleration voltage), only times the mean value 
of the z-component of operators will remain. However, 
for example, in 2p — > Sd transitions of iron, cobalt or 
nickel at 200 kV acceleration voltage the z-component 
of q- vectors is approximately 1/3 of the reciprocal lat- 
tice vectors, i.e. non-negligible. On the other hand, in 
systems with weak spin-orbit interaction or in strong 
enough external magnetic field, the magnetic moments 
in the sample will have a direction dictated by the ex- 
ternal magnetic field, which is usually the optical axis of 
the transmission electron microscope (the z axis). Then 
again one can approximate Q • (L) « Qz{Lz)- 

If none of these assumptions hold, the sum-rules give 
a projection of the spin, orbital and magnetic dipole mo- 
ments on the direction of Q. The following reasoning 
then helps to interpret the sum- rules predictions. The 
most common EMCD geometries can be very well simu- 
lated using the systematic row approximationjii There- 
fore all q vectors lie within one plane in the reciprocal 
space and a vector product of any two of them is normal 
to this plane - i.e. all Q vectors have the same direction. 

An important feature of the EMCD sum-rules is the 
fact that these allow in principle any detector position 
to measure magnetic moments. I.e. the experiment is by 
no means limited to placing the detector on the Thales 
circle above two strongest Bragg spots.— This observa- 
tion is particularly important for the new experimental 
geometries^ and for the optimization of the signal to 
noise ratioi^^ 

The accuracy of the derived sum-rules naturally needs 
to be evaluated. We have investigated this by com- 
paring the spin and orbital moments obtained from 
the EMCD sum-rules for a complete simulation of the 
dichroic signalii of bee Fe, hep Co, and fee Ni, to ex- 
act spin and orbital moments. The latter values are 
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FIG. 1: (color online) Comparison of ab initio orbital (blue 
solid line), spin and magnetic dipole moments (red dashed 
line) to the simulated sum-rules prediction (symbols) as a 
function of effective U parameter (see text). 

calculated from the quantum mechanical expectation 
value of the spin and angular momentum operator, us- 
ing wavefunctions from first principles theory. In the 
past such calculations have been shown to agree very 
well with measured moments of e.g. Fe, Co, and NiJ^ 
In our calculationaii we have employed the strength of 



the electron-electron interaction (represented by the so- 
called Hubbard U parameter) as a means to vary the size 
of the spin and orbital moments. Figure [1] demonstrates 
an excellent agreement between the moments extracted 
from the sum-rules and the calculated operator mean val- 
ues. It should be noted that in general the value of (T^) 
is low, but for hep Co in the range 0.2 < U < 0.4 Ry it 
attains a rather large value. This is caused by changes 
in the electronic structure and will be analyzed in detail 
elsewhere. 

Finally we note that an actual application of the 
EMCD sum-rules to extract absolute number of spin and 
orbital moments is more difficult than in XMCD exper- 
iments. Due to the effects of the dynamical electron 
diffraction the sum-rule analysis relies on knowledge of 
the sample thickness dependent function, which enters 
as a different prefactor within imaginary and real parts 
of the MDFF's, respectively. This can be calculated us- 
ing the dynamical diffraction theoryJ^ 

On the other hand, the ratio of spin and orbital mo- 
ments can be extracted in a completely analogous way as 
in the XMCD, because both spin and orbital momentum 
sum-rules depend on the same prefactor (up to a factor 
3/2). Similarly if the 'extra terms' in N and spin-orbit 
sum-rules are negligible, the ratio of (L • S) and Nf, can 
be obtained. 

In conclusion, we have derived a complete set of sum- 
rules for the general ELNES experiment. Our theory 
provides a foundation for the use of the EMCD tech- 
nique as a quantitative tool in magnetic studies. Such 
quantitative EMCD analysis will open up new avenues 
for research in the field of nano-magnetism, owing to an 
excellent spatial resolution in combination with atom- 
specificity, primarily for obtaining information about 
spin and orbital moments. 
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Electronic structure calculations were done using the FP- 
LAPW method with 40000 or 20000 k-points in the first 
Brillouin zone for Fe, Ni or Co, respectively. Atomic sphere 
radii were set to 2.2 Bohr radii for Fe and Ni and 2.3 Bohr 
radii for Co, respectively. Approximately 100 basis func- 
tions per atom were used. 



